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Abstract 



Classical random matrix ensembles with orthogonal symmetry have the property that the joint 
distribution of every second eigenvalue is equal to that of a classical random matrix ensemble 
with symplectic symmetry. These results are shown to be the case r = 1 of a family of inter- 
relations between eigenvalue probability density functions for generalizations of the classical 
random matrix ensembles referred to as /^-ensembles. The inter-relations give that the joint 
distribution of every (r + l)-st eigenvalue in certain /3-ensembles with (3 = 2/(r + 1) is equal 
to that of another /^-ensemble with /3 = 2(r + 1). The proof requires generalizing a conditional 
probability density function due to Dixon and Anderson. 
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1 Introduction 



1.1 The setting and summary of results 

The Dixon-Anderson conditional probability density function (PDF) refers to the function of 
{Xj} specified by [3l, T] 

r(.i) • • • r(.„) nr,,<..n(«. - «^)^^"^-^ M M ' ^ " ^' ^ ^ ^ 

where generally = 1 if condition ^ is true, XA = otherwise, and here condition A is the 
inequalities 

ai > Ai > 02 > A2 > • • • > A„_i > an 

specifying an interlaced region. An analogous conditional PDF for angles {tpj}j=i,...,m due to 
Forrester and Rains [13], is specified by 

where On = Itx and R denotes the interlaced region 

< V'l < 6*1 < V2 < 6*2 < • • • < V'n < 27r (1.3) 

(in fact in [13] the parameters were specialized to ai = • • • = On-i = a, but the working 
therein applies without this specialization). The conditional PDF (jl.ip is fundamental in the 
theory of the Selberg integral [91 Ch. 3], [TJ]. It permits various generalizations, which are 
intimately related to Macdonald polynomial theory [T7j and the theory of multivariable elliptic 
hypergeometric functions [18]. 

It is the purpose of this paper to introduce different generalizations of (II. ip and ()1.2p . and 
to discuss their consequences to random matrix theory. To state our generalizations, which will 
be derived in Section 2, we require for the normalizations the Selberg integral [20j 



5'Ar(Ai, A2, A) := 



/ dh--- dtNW *^ (1 - *')^' n 1*^ - *j 

•'^ -'^ 1=1 l<j<k<N 



N 

|2A 
(.J I 

l<j< 

"^^r' r(Ai + 1 + jA)r(A2 + 1 + jA)r(l + (j + 1)A) 
M r(Ai + A2 + 2 + (iV + j-l)A)r(l + A) ' 



and the Morris type integral (see [9l Ch. 3]) 

-1/2 1-1/2 N 



MN{a,X) := (27r)^ / d^i • • • / d^iv H 1 1 + e^^f " J] {e'^" - e 

■^-^1'^ J-'^/'i 1=1 l<j<k<N 

-- Aj + 2a + l)r(A(j + l) + 
(r(Ai + a + l))2r(l + A) 



idj |2A 



sTV TT r(Aj + 2a + l)r(A(j + l) + l) 



Theorem 1. Let r G Z+. T/ie Dixon- Anderson PDF U.l\) is the r = 1 case of the family of 
conditional PDFs 

1 Ui<j<k<r(n^i)i^j - ^k?^^"^^^ 'n'^TTlA a\^p-K, (16) 

T^TT 7 \r(s-+si-2/(r+l)) 11 lll'^J""?! U-Oj 
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where Ar is the interlaced region 

aj > K{j-^i)+i > K(j-i)+2 > ■ ■ ■ > K{j-i)+r-i > flj+i (j = 1, • • • , - !)• (1-7) 
and the normalization C is specified by 

n— 1 ^ I 

^=Il^_^r{^Sp + 2{l- l)r/(r + l)-l,si- 1, l/(r + 1)) . (1.8) 
1=1 ' p=i 

The circular analogue of the Dixon- Anderson PDF is the case r = 1 of the family of 

conditional PDFs 

- ^^i<j<k<rn\^ ^ TTTTle^'^^ -e^^H^^-^R (19) 

"-^ Hl<j<fc<n P ^ I j=lp=l 

where = 27r, i?,. denotes the interlaced region 

6j_i < '0(r-l)i+l < V'(r-l)i+2 < • • • < V'rj < % (j = 1, • • • , ?^) (1-10) 

wii/i ^0 0, and the normalization C is specified by 

C = (7|{,^}^|„^|M, ( 2 ( E + - + 1) - ^) , + 1)) • (1-11) 

p=i 

As to be revised in Section 3, both (jl.ip and (jl.2p have consequences in random matrix theory, 
in that they provide inter-relations between ensembles with orthogonal, unitary and symplectic 
symmetry. Likewise, it will be shown in Section 4 that for general r G Z"*" the generalizations 
(|1.6|) and (|1.9p have consequences in random matrix theory. To state these consequences, let 
ME^^jv((7) refer to the matrix ensemble specified by the eigenvalue PDF 

1 ^ 

cUsi^i^ n \^k-x,f (1.12) 

1=1 l<j<k<N 

(unless otherwise stated, C will denote some normalization), and let the eigenvalues in (jl.l2p 
be ordered 

Xi > X2 > • • • > xn- (1-13) 

Also, let CE^ ^ refer to the so called circular Jacobi /3-ensemble of unitary random matrices 
(see e.g. [Hi Ch. 2]), which is specified by the eigenvalue PDF 

lfj|l-e^^f n |e'''=-e*'^|^ (1.14) 

1=1 l<j<k<N 

and suppose the angles are ordered 

< 01 < 02 < ••• < ^iv < 27r. (1.15) 
Denote by Dp(ME^ 7v(5')) the joint margiiicLl distribution of x^, X2p5 ^Sp? • • • in 

(fLT2D . and de- 
note by Dp(CE^^) the joint marginal distribution of 0p,92p,0^p, . . . in (11.141) . Following the 
nomenclature of |12| . the D here stands for a decimation procedure. 
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Theorem 2. Let r G 



One has the inter-relations between matrix ensembles 



Dr+l{ME2/(r+l),ir+l)N+rix"'{l - xf) 

D,+i(ME2/(,+i),(,+i)jv((l-a:)^) 

Dr+l (ME2/(r+l),{r+l)N+r (x^e""") 
Dr+1 (ME2/(r+l),(r+l)Ar (e""") 
Dr+l(ME2/(r+l),(r+l)Ar+r(e~^ ) 

Dr+l iG'i^2/{r+l),{r+l)N+r 
Dr+1 (CE2/(j,_,_l)^(,r_,_l)^ 



ME. 



2(r+l 



(r+l)a+2r. 



1 



N(r+l)fe+2rN 



ME2(,+i),jv((l-x) 



(r+l)fe+2r\ 



ME; 

ME. 



ME. 
CE 



2(r+l) 

2(r+l) 

2(r+l) 
(r+l)b+2r 
'2(r+l),Af 



,N(.e 
,N{e 



{r+l)a+2r ^-{r+l)x-^ 
-{r+l)x\ 



-{r+l)x2 



(1.16) 



In Section 5 these inter-relations, applying at the level of joint eigenvalue PDFs, will be shown 
to imply analogous inter-relations linking /? = 2{r+l) and /? = 2/(r+l) for marginal distributions 
of single eigenvalues. Depending on these single eigenvalues being in the neighbourhood of the 
spectrum edge, or in the bulk, there are three distinct — > oo scalings, implying corresponding 
inter-relations for marginal distributions in the scaled ensembles. 



1.2 Significance of the P ensembles 



One interpretation of (jl.l2p is as the Boltzmann factor for a classical gas at inverse temperature 
P with potential energy 



1 ^ 



^<j<k<N 



Because of the pairwise logarithmic replusion (two-dimensional Coulomb law), such a classical 
gas is referred to as a log-gas. This interpretation allows for a number of properties of correlations 
and distributions to be anticipated using arguments based on macroscopic electrostatics. Of 
interest for purposes of the present study is a prediction relating to the probability E^"^^{n; 2t) 
that an interval of size 2t in the bulk contains exactly n eigenvalues, in the N ^ oo limit of 
(jl.l2p , scaled so the eigenvalue density is unity (by universality, for a large class of g this limit 
is expected to be independent of g). The prediction is that [HE] 



logE^0-'Hn;2t) 



t.n — >cxD 



l]^t+^l-^-^ 

2\ 2 2 



log + 1 

n 



(1.17) 



We will see in Section 15.21 that this is consistent with the final relation in p.l6|) . 

A second interpretation of (jl.l2|) is as the absolute value squared of the ground state wave 
function for certain quantum many body problems with inverse square pair potentials (Calogero- 
Sutherland systems). For example, with g{x) = e~^^ (jl.l2p is equal to iV'oP where iJjq is the 
eigenfunction corresponding to the smallest eigenvalue of the Schrodinger operator 



N 



32 N 



H 



j=l ^^j ^ j=l l<j<k<N 



Such interpretations are restricted to a total of four examples of ()1.12p or (jl.l4p (see e.; 
namely g{x) equal to 



mi 



x^il-xf (0<x<l), x'^e-'^ (x>0), 
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in (|1.12p and 6 = in ()1.14p . In the study of these ground states, identities of a different type 
to Theorem [2] relating /3 to 4//3 (/? even) have previously been encountered. These are so called 
duality relations, an example being [2] 



N 



IK* 



i=i 



IX 



ME, 



'4//3,m 



The four examples of (jl.l2p or (I1.14p which permit interpretations as the ground state of 
quantum many body systems can also be realized as the eigenvalue PDF for certain ensembles 
of random matrices. Thus, for example, with denoting value drawn from the square root 
of the gamma distribution r[A;/2,l], and N[0, 1] denoting a number drawn from the standard 
normal distribution, the tridiagonal matrix 



f/3 : = 



N[0,1] X{N-i)f3 
X(iV-i)/3 N[0, 1] 

X{N-2)I3 



N[0, 1] X(N-3)f3 



X2/3 



N[0,1] 



Xl3 

N[0,1] 



has its eigenvalue PDF given by ME^^Tvle-'^'/^) g]. Taking various scaled N ^ oo limits 
(see Section 15. 2p of these ensembles of random matrices leads to a description of the limiting 
eigenvalue distributions in terms of stochastic differential operators [3 \T9\ I16j . Further, for 
eigenvalues in the bulk, there is a description in terms of a process associated with Brownian 
motion in the hyperbolic plane |21j . 

It is an open problem, for general r, to derive the results of Theorem [2] as consequences of 
the eigenvalue PDFs being realizable from concrete random matrix ensembles. In the case r = 1 
such matrix theoretic derivations have been given (excluding the sixth identity which relates to 
CE^) using matrix realizations not applicable for general r \12\ I13j. 



2 Proof of Theorem [T] 

2.1 Analytic continuation of L^„({ap}) — a cancellation effect 

Let 



L. 



■,ni{ap}) := dXi ■ ■ ■ dXr(^n-i) n (-^J 

l<j<k<r(n-l) 



r{n~l) n 

j=i p=i 

(2.2) 



^r,n({ap}):= n (a,-a,)^(^^+^'=-2/('-+i)), 

l<jr<A:<n 

where in (j2.ip Ar is given by (jl.7p . To prove that (II. 6p is a PDF in {Xj} we must show that (12. ip 
is proportional to (j2.2p . In preparation for this task, we will study the analytic continuation 
of (j2.ip . and in particular its value upon the interchange of the order of {op}. For the latter 
purpose, it suffices to consider the effect of the interchange of an arbitrary neighbouring pair of 
parameters a/, a^+i say. One finds that, up to a phase, the analytic continuation is a symmetric 
function in {{sp,ap)}. 
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Proposition 1. Interchanging the order of the {flp} on the real line via analytic continuation, 
and making a corresponding change to the order of {sp}, gives back the same integral represen- 
tation \2. up to a phase. In particular 

LrAWp})\^^^.^^^ = e— ^(^'+^'+-2/(^^+i))L,,„({a,}). (2.3) 

The proof of Proposition [T] relies crucially on a certain cancellation effect, isolated by the 
following result. 

Lemma 1. Consider an arrangement of r O's and q I's {1 < q < r) in a line. Let this be 
considered as the sequence A = {nj)j=i^...^r+q with each Uj =0 or 1. Further let 

K{n) = [^' ""^^^ (2 4) 

^Vh) \ #o'stotherightofnj, = 1, ^ > 

and use this to specify the statistic 

r+q 

K{A) = Y,K{n,). (2.5) 



One has 



Y^^-2.^K{A)/{r+l) (2.6) 



A 

Proof. The definition (12.41) can be written 



r+q 

K{nj) = Uj (l-rij), 
k=j+i 



and this substituted in (j2.5p gives 



r+q r+q 

^(^) = XI ^^^^ + Q-j)-^ ^j^k- 
j=l j<k 

Noting 

r+q ^ r+q ^ ^ r+q ^ 

j<k j=l j=l 

we see that ()2.6p is equivalent to showing 



^g-27riE;i?nj(r-+g-j)/(r+l) (^2.7) 
A 

Now, the LHS of (|2.7p is equal to the coefficient of z'' in 

F{z):= Yl e-2-^;;?-i(^'-+''-^')/(^'-+i)z^.-i"^ (2.8) 

ni,...,nr + q = 0,l 

The multiple sum factorizes into a product of r + g single sums, giving the factorization formula 
F{z) = llil + ^e-2-(j-i)/(-+i)) = (1 _ + ^e-2.riO-l)/(r+l)^^ 
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where in obtaining the second equahty use has been made of the simple identity Hz^iCl 
2g27r«(/-i)/Af-j _ ^ _ ^ "Yy-g read off from this that for some ci, . . . , Cq-i, 



F{z) = 1 + ciz + • • • + Cg^iz'^-' + (-!)'■ z^+i + (-l)^ciz^+2 + . . . + {-lycg.iz'+'i. 

In particular (recalling that 1 < q < r) there are no terms proportional to z"^, and so (j2.7p must 
hold. □ 

2.2 Proof of Proposition [1] 

As written (|2.ip is only defined for real {a^}, ordered so that ai > • • • > a^. It can be defined 
for general complex numbers by analytic continuation. This requires first replacing the absolute 
value in (12. ip according to 



Xj — ap, for Xj > Op 
ap — Xj , for flp > Xj , 



so that the integrand consists entirely of power functions. According to Cauchy's theorem, 
the contours of integration can now be deformed into the complex plane, provided no contour 
crosses a branch cut of the power functions. With the latter specified by z°' := \z\°'e^°'^^^^ , 
—TT < argz < TT, the branch cut is on the negative real axis of the complex z-plane. Moving 
the position of the Op's into the complex plane in this setting gives the analytic continuation of 

Our interest is in this analytic continuation when ai and ai-^i swap places on the real axis 
for general I = 1,2, ...,n — 1. This analytic continuation corresponds to the limit that the 
contours tend to the real line in the second configuration of Figure [H and the endpoints {dj} 
are appropriately related to the endpoints {aj}. 

m m m m 1631 clXlS 



I real axis 

ai-i 

Figure 1: The contours from a;+2 to a;+i, a;+i to ai, and ai to a;_i are deformed to the contours 
joining the corresponding tilded variables. Our interest is in the limit that dj = aj {j ^1,1 + 1), 
di = a/+i, di^i = ai and all contours in the second diagram run along the real axis. In the case 
I = n — 1 the contour from dij^2 to a/+i is to be deleted, while in the case / = 1 the contour from 
di to di-i is to be deleted. 

To study the integrand of (|2.ip in the case of the second configuration, for notational con- 
venience set A(j_i)j,_|_i, = A^'^^ (z/ = 1, . . . , r), which are to be referred to as species j. Then for 
the first configuration 

a,+i < Xf < Af < • • • < Af <a, (j = 1, . . . , n - 1), 



^1+2 ^1 
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while the integrand itself reads 



j=ll<u<fj,<r l<j<k<niy,fi=l 



r n n 

X 

iy=l j=l p=l 



r n ri 

( n n n i^r - =■■ ^rAi^f^h {ap}). (2.^ 



In the second configuration let the integration variables of (I2.9P be tilded so that A^'''' i— > X^j^^ . 

On each contour from dj to Oj+i the integration variables {X^^^}i^=i,...,r are ordered so that 

A^^^ >- X^^^ >-■■■>- Xj~\ where the ordering >- is the descending order induced by the direction 
of the contour. In the limit that all the contours run along the real axis, this ordering corresponds 
to 

dj > Af ) > X^ >■■■> X^^ > j = 1, . . . ,n - 1 (j / 0, 

di+i > Af'^ > Af'~^^ > • • • > a[^^ > di. (2.10) 

Furthermore, we see from Figure [1] that for some p = 0, . . . , r, g = 0, . . . , r, 

ai+i > a[|\ > Xl% >■■■ aJJ, > di (2.11) 

a,+i > a;:7^+^) > a;:7+') > • • • > x^], > a, (2.12) 

(note that these configurations are empty if p = 0, q = respectively). In other words, between 
di and 0/+1 there are r coordinates of species I, p oi species I + 1 and q of species I — 1. 

A crucial feature of the contour integrals is that only configurations with p = q = con- 
tribute, due to cancellation effects for p and/or q non-zero. To quantify the latter, consider first 
the case that p = while g > 1, and suppose that to begin the r species I variables are to the 
left of the q species I — 1 variables in the interval (a;, a;+i). We see from (|2.9|) that interchanging 
the position of coordinates corresponding to different species does not change the magnitude of 
the integrand but it does change the phase, with each interchange of a species / — 1 and left 
neighbouring species I contributing _ Hence for a general ordering of the r species / 

variables and q species I — 1 variables amongst a given set of (r + q) positions in (a;,a;+i) the 
phase is given by 

g-27riX(.A)/(r+l)^ 

Here K{A) is as in Lemma[T]with the O's corresponding to species / and the I's to species I — 1. 
But Lemma [J tells us that if we sum ()2.13p over all arrangements we get zero, which is the 
claimed cancellation effect in this case. 

Essentially the same argument, making use of Lemma [1] with the role of the O's and I's 
interchanged, gives cancellation of the contribution to the contour integrals from configurations 
with g = 0, p > 1. It remains to consider the cases p,q > 1. In such cases, with the positions 
of the species / — 1 coordinates fixed (we could just as well fix the position of the species I + 1 
coordinates), we see that the contribution to the phase of each such coordinate is equal to 
g-27nZ*/(r+i)^ where I* is the number of both species I, / + 1 to its left, and in particular is 
independent of their ordering. But we know that summing over this latter ordering gives the 
cancellation (j2.ip . so in all cases there is no contribution from non-empty configurations (|2.1ip 
and dsn. 
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As a consequence of both ()2.1ip and (j2.12p having to be empty for a non-zero contribution 
to the contour integral, it follows that (j2.10p can be supplemented by the requirements that 

ai > > > ■■■ > a«+2 
di-i>~x['\>~x[^\>---\\'\>di+i. 

Up to a phase, this contour integral is precisely (|2.ip with the position of ai and a^+i in- 
terchanged, and correspondingly si and s^+i interchanged. The phase is straightforward to 
calculate, giving as the final result (j2.3p . □ 



2.3 Proof of Theorem 1 for (fTTej) 

As remarked below (|2.2p . we must show that Lr,n({ap}) is proportional to Rr,ni{0'p}), and then 
determine the proportionality. For the former task, our strategy is to show that Lr^n{{cLp}) 
factorizes into a term singular in {cp}, and a term analytic in {a^}. The singular factor is 
precisely Rr,n{{0'p}) , while a scaling argument shows that the analytic factor must be a constant. 
Intermediate working relating to the singular terms allows the proportionality to be determined. 

Consider I/r,n({ap}) as an analytic function of ai in the appropriately cut complex ai-plane. 
Singularites occur as ai approach any of a2, . . . , a„. The singular behaviour as ai approaches 02 
can be determined directly from (|2.ip . Thus, as ai — > a2 the integral over species 1 effectively 
factorizes from the integral over the other species, showing 

n 
p=3 

xLr,n-l{{ap}p=2,...,n)\si^si+S2+2/{r+l)-lP{C'l " {ap}p=2,...,n) (2-14) 

where F(z; {ap}p=2,...,n) is analytic about z = and equal to unity at z = 0, and 
Iriai,a2) := / dX^^^ ■ ■ ■ dX^{^ 

X n (AS"^-Al^^)'/("+^)n(ai-Ai")r-i(A;'^^-a2r-^ (2.15) 

l<!/</x<r i/=l 

Thus the singular behaviour is determined by the singular behaviour of 1^.(01,02). This in turn 
is revealed by a simple scaling of the integrand, which shows 

Ir{aua2) = (ai - a2r^''-'^/^'+''>+'-^''+''~''> ^Sr{si - 1, ^2, l/(r + 1)) (2.16) 

where ^^(Ai, A2, A) denotes the Selberg integral ()1.4p . 

For the singular behaviour as oi approaches (A; 7^ 2), we make use of Proposition [1] 
which says that up to a phase the function of {ap}p=i^...^„ obtained from Lr,n({op}) by analytic 
continuation is symmetric in {(flp, Sp)}. Hence as a function of oi it must be that 



n 



LrAWp}) = lliaj - afc)'^('-i)/(^+i)+^(^^+^'=-i)F(ai; {ap}p=2,...,n), (2.17) 



k=2 



where F is analytic in ai. Further, repeating the argument with i^r,ra({flp}) regarded as a 
function of 02 , . . . , fln in turn shows 

UAiap}) = RrA{ap})G{{ap}) (2.18) 
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where G is analytic in {op} and symmetric in {(a^, Sp)}. 

It remains to determine G. This can be done by considering the scaUng properties of both 
sides of (|2.18p upon the replacements {op} ^ {ccp}, c > 0. After changing variables Afc i-^ c\k 
(A; = 1, . . . , n(r — 1)) in (j2.ip we see 

L.,„({cap}) = cr(n-l)+r(n-l)(r(n-l)-l)/(r+l)+r.(n-l) E^^iC^p-D^^^^d^^}) (2.I9) 

while we read off from (j2.2p that 

iir,n({cap}) = ( n c^(^^+^'=-'/('-+^)))/?.,„({ap}). (2.20) 

l<j<fc<n 

Straightforward simplification shows that the exponents of c in both equations are in fact equal, 
and so Lr,n({ap}) and Rr,n{{ap}) are homogeneous of the same degree. Thus according to ()2.18p 
the function G{{ap}) must be homogeneous of degree 0. Because G is analytic in {flp}, this 
requirement implies that G is actually independent of {op}, so Lr,n({ap}) is proportional to 

Rr,n{{a-p})- 

The remaining task is to compute the proportionality. For this purpose, note that because 
the normalization is independent of the Op's, 

lim Lr,n{{aj})/Rr,n{{aj}) = C. (2.21) 
ai ,...,an — 

On the other hand, iteration of ()2.14p with the substitution (I2.16P allows this limit to be com- 
puted in terms of Selberg integrals, giving the result (jl.Sp . □ 



2.4 Outline of the proof of Theorem [U for ( fTTOl ) 

The strategy used to establish Theorem [T] in the case of (11. 6p requires only minor adjustment 
to also establish Theorem [1] in the case of ()1.9p . Thus with 

/rn n 
#i---rfVrn n |e*'^'= -e^^^f/(^+^) J] JJle^'''^ -e^^^-r^-i (2.22) 
'^Rr l<j<k<rn j=lp=l 

and Wp := e^^^, the analytic continuation in {wp} is constructed by rewriting the absolute values 
in the integrand according to |n — = {u — v){l/u— 1/v), writing zj = e*^J {j = 1, . . . , rn)) and 
considering (j2.22p as a contour integral over arcs of the unit circle in the appropriate complex Zj- 
planes. By then appealing to Lemma[T]the analogue of Proposition [T] can be proved, establishing 
that the analytic continuation of L{{'Wp}) is a symmetric function of {(tfp, ap)}p=i,...,n- 
The next task is to establish the factorization 

QrAi^p}) = SrA{^p})G{{Wp}) (2.23) 

where 

^ /r IN TT // Wi\/' till.. W '•((aj+afc)/2-l/(r+l)) 

i<j<k<n •' 

while G{{'Wp}) is a symmetric function of {{wp,ap)} analytic in {wp}. This is done by first 
studying (j2.22p in the limit wi — > W2 to obtain the analogue of ()2.17p . then using the analogue 
of Proposition [1] to deduce (|2.23p . But both (5r,n and Sr,n are homogeneous of degree zero, and 
so G must therefore be a constant. 
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The constant C is independent of the OiS and so can be computed according to 

hm QrAW'^))/SrAW'^))=C. (2.25) 

Ol,...fin-\^0 

Now it fohows from (p:22]) that for 6*1 ^ (92, 

n 
p=3 



a2i-^ai+02+2r/(r+l) — 1 



where Ir is specified by (j2.15p and evaluated in terms of the Selberg integral by ()2.16p . Iterating 
this and recalling On = 2tt allows the limit (j2.25p to be computed in terms of products of Selberg 
integrals and a Morris integral, giving (jl.lip . □ 

3 Random matrix consequences of the case r = 1 

It has been pointed out in [12^ Section 4.1] that special cases of (jl.ip permit random matrix 
interpretations relating to the Jacobi random matrix ensemble, specified by the eigenvalue PDF 

1 ^ 

-Jlx1{l-xif W \xk-x/, 0<Xi<l, (3.1) 

1=1 l<j<k<N 

or equivalently, with the ordering (jl.lSp assumed and it implicit that the support of x°(l — x)^ 
is [0,1], MEp^n{x°'{^ — x)^), in the cases /3 = 1,2 or 4. These cases can be realized in terms 
of certain random matrices with Gaussian entries (see e.g. [9]), the entries being real, complex 
or quaternion real for /? = 1,2 or 4 respectively. The corresponding PDF on the matrices 
is then invariant under similarity transformation by orthogonal (/? = 1), unitary (/? = 2), 
unitary symplectic (/? = 4) matrices, giving rise to the alternative notations MEi^tv ^ OEjv, 
ME2,Ar 1-^ UEat, ME4^Ar i— > SE^r, which are to be used below. Following [IT], let us also change 
the name of the special decimation procedure D2, calling it instead by the word even. In terms 
of these notations, it has been proved in [11] that 



evenOE27v+i(a;^""^^/^(l-x)(''-^)/2)) = SENix''+\l - xf^^) (3.2) 

evenOE27v((l-2;)(^~^)/2) = SEw((l - x)^+i). (3.3) 

Consider first the derivation of (j3.2p from (jl.ip . For this set 

n = N + 2, ai = 1, ajv+2 = 0, si = (6+l)/2, sn+2 = (a + l)/2, sj = 2 {j = 2,...,N + l). (3.4) 

Noting that then 

n—l n—1 n 

iiaf-'^/\i-a,f-'y^ n («.--«^) n (A.-A.)nni^.--«pi^^"' 

j=2 l<j<k<n l<j<k<n-l j=l p=l 

ocOE2;v+i(x("-^)/2(l-x)(^-i)/2)) (3.5) 



while 



ra-l 

JJa('^-3)/2(i_a^.)(fe-3)/2 Yl (a,--afc)^^+^^ocSE^(x'^+i(l-x)^+i) (3.6) 

j=2 l<j<k<n 
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we see that ()3.2p follows from (jl.ip by integrating over {Xj} and using the fact that with respect 
to the latter variables (jl.ip is a PDF and so integrates to unity. To derive (j3.3p from (jl.ip . set 



n = TV + 1, ai = 1, si = (6 + l)/2, = 2 (j = 2, . . . , A^). (3.7) 

We then have 



-1 



71 71 — 1 7i 

n(l-a,)(^-3)/2 J] (a, -a,) J] " ^'c) J] 11 l^^' 

i=2 l<j<fc<n l<j<A:<n-l j=l p=l 



while 



na-a,)^'"'^/' n (a,-afc)^^+^'=ocSE^((l-x)^+i) (3.9) 

j=2 i<j<k<n 

and so p.3p is also seen as a consequence of (jl.ip being a PDF and so integrating to unity. 

The inter-relations (j3.2p and (j3.3p have companions, which in ^llj were proved in fact to be 
equivalent to the originals. These state that 

even(0Eiv(x("~i)/2(i _ ^)(fe-i)/2) y 0E^+i(x("-i)/2(i _ ^)(fc-i)/2))^ = UEjv(x'*(l - x^)) (3.10) 

even^OEjv((l — x) ("-1)/^) U 0E^(1 - x)(''-i)/2)) = UE^((1 - x')). (3.11) 

In general, for two matrix ensembles ME^^^„ and ME^j^m; the operation ME^^_„UME^2,m denotes 
the ensemble of (n + m) eigenvalues formed by sampling independently from ME^j^„, ME^j,™, 
superimposing the resulting eigenvalue sequences, and labelling from right to left. We know 
from [llj that the eigenvalue PDF of OE„(/)uOE„-(_i(/) is proportional to 

2n+l 

-1 - X2k-l) Y\. (^^2j - X2k) (3.12) 
1=1 l<j<k<n+l l<j<k<n 

while the PDF of OE„(/)uOE„(/) is proportional to 

2n 

-1 - X2k-l)ix2j - X2k)- (3.13) 

1=1 l<j<k<n 

Given these facts, the inter-relations (j3.10p and (j3.1ip can be understood as corollaries 
of (jl.ip . Explicitly, to deduce (j3.10p . choose the parameters as in (j3.4p except that Sj = 1, 
(j = 1, . . . , N). Then, recalling (j3.12p . we have that the LHS of (j3.5p is proportional to the LHS 
of (1330]) . while the LHS of dM]) is proportional to the RHS of (iSTOl) . Similarly, to deduce (IHTTT) 
from (jl.ip . choose the parameters as in (j3.7p but with Sj = 1 (j = 2, . . . , N). Recalling (j3.13p . 
we then have that the LHS of (j3.8p is proportional to the LHS of (j3.1ip . while the LHS of (j3.9p 
is proportional to the RHS of (j3.1ip . 

The Jacobi /3-ensemble (j3.ip permits two well known limiting cases. Thus, with the replace- 
ments Xj Xj/b and taking the limit 6 — > oo, one obtains the eigenvalue PDF for the Laguerre 
/3-ensemble, 

1 ^ 

-J]<e-^'' n l^'k-xjf, 0<xi<oo. (3.14) 

1=1 l<j<k<N 
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Also, with a = b and xj ^ ^(1 — Xj/^/2B), taking the hmit 6 — > oo one obtains the eigenvalue 
PDF for the Gaussian /^-ensemble 

^Yi^-^'/' n 1-^-^/- (3-15) 

1=1 l<j<k<N 

The inter-relations (j3.2|) . and (|3.1U|) . permit these same limiting cases. For the 

Laguerre limit we read off from these that 

evenOE27v+i(2;("-^)/^e-^/2) ^ SEN{x^+^e-^) 
evenOE2jv(e~''/^) = SEtvIc"^) 

even (OE7v(x('^-^)/2g-x/2) ^ OEjv+i(x('^-i)/2g-x/2^^ ^ UE^(x'^e-^) 

even (0E^(e-^/2) ^ OEiv+i(e-^/^)) = UE^(e-^), (3.16) 

while only (j3.2p and ()3.10p permit Gaussian limits, which read 

evenOE27v+i(e-"'/2) = SEiv(e-"') 

even (0E7v(e-^'/^) U OE;v+i(e"'''/^)) = UE^(e-^'). (3.17) 

These were first obtained in [llj in the context of classification theorems relating to the form 
of f{x) in 0E7v(/) which permits such identities. We remark too that as observed in [12|, the 
Dixon-Anderson PDF (jl.ip permits analogous Laguerre and Gaussian limits, and these have 
(I3.16P and ()3.17p as corollaries. 

We now turn our attention to the revision of some consequences in random matrix theory 
relating to the conditional PDF ()1.2p [13j. For this we recall that the circular ensembles COE„, 
CUE„, CSE„ of unitary random matrices have their eigenvalue PDF proportional to 

Y[ |e*^*-e^^^f (3.18) 

l<j<k<n 

for P = 1, 2, 4 respectively. Further, we require knowledge of the fact that the PDF of COE„UCOE„ 
is proportional to [15] 

II sin(02fc - e2j)/2 sin(02fe-i - 02j-i)/2 (3.19) 

l<j<A:<n 

with {62j} and {6*2^-1} interlaced according to 02j-i < d2j < ^2j+i (j = Ij • • • j J^)) ^2n+i := 27r. 
In ()1.2p regard 9 variable with 

ipn < On, 6l„mod27r < Vi- (3.20) 



Let R denote the interlaced region (jl.3p supplemented by (j3.20p . Then, from the fact that (|1.2 
is a PDF in n we have 



ii>l \aj-l 



„ n n 

•'^ l<j<k<n j=il=l 

oc Yl |e*^^ -e^^^ |"^+°^ (3.21) 

l<j<k<n 
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Setting aj =2 (j = 1, . . . ,n) we obtain the inter-relation between circular ensembles [6] 

altCOEsn = CSE„, (3.22) 

where the operation alt refers to the distribution of every alternate (second) eigenvalue as ordered 
around the circle. Further, after making use of (j3.19p . setting aj = 1 (j = 1, . . . , n) we obtain 
the companion inter-relation 

alt (COE„ U COE„) = CUE„. (3.23) 

An ensemble which generalizes (j3.18p is the circular Jacobi /3-ensemble (|1.14p . which in the 
special cases of /3 = 1,2,4 is to be referred to as COE^, CUE^, CSE^ respectively. It follows 
from (jl.2p itself (and thus e*^" therein equal to unity) that 

evenCOE^^+i = CSE^^+2, even(COE^+i U COE^) = CUE^^+\ (3.24) 

The first of these contains ()3.22p as a special case. To see this, note from ()1.14p that in general 

C^ftAf °^ ^^/3,Af+llfiV+i=0- (3.25) 

Setting 6 = 1 in the first formula of ()3.24p is then seen to reclaim (I3.22p with n = N + 1. 

4 Random matrix consequences for general r 

In this section Theorem [2] will be established. Consider the integrand (j2.ip in the case 

n = N+2, ai = 1, aN+2 = 0, si = 6+1, sn+2 = a+l, sj = l+2/(r+l) (j = 2, . . . , iV+1) (4.1) 

(cf. (j3.4p ). After multiplication by a suitable function of {op} we see that a matrix ensemble 
PDF of the form (jl.l2p results. Explicitly one has 

J|a-2/('-+l)(i_„^y-2/(.+l) -Q _ „^)2/(.+l) -Q _ ^^)2/(.+l) 

i=2 i<j<k<n l<j<k<r{n-l) 

r{n-l) n 

X n ni^J-'^pl"'"^ ^^^2/{r+lUr+l)N+r{x''{l-xf). (4.2) 
j=l p=l 

Multiplication of (j2.2p by this same factor gives 

n-l 

n a;-'/(''+^) (1 - a,f-'/(''+^) n - afc)'^(«^+^^)-2(--i)/('-+i) 

j=2 ^<j<k<n 

oc ME2(,+i),jv(2;("+^)"+^"(l - x) ('■+1)^+2'-). (4.3) 

The fact, as implied by Theorem [1] as it relates to ()1.6p . that (j2.ip and (j2.2p are proportional, 
tells us that integrating (j4.2p over the variables {Aj}, (14. 3p results. Now, by the ordering 
implied by (|1.7p . and the fact that ai is fixed, integrating over the variables {Aj} gives the joint 
marginal distribution of the variables in positions (r + 1), 2(r + 1), 3(r + 1), . . . as counted in 
a descending order. Thus the integration in (|2.ip is the operation D^+i, and so the equality 
between integrating ()3.2p over {Aj} and (13. 3p can be expressed as the first of the relations in 
(fTTHl) . 
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To derive the second relation, set 

n = N + 1, ai = I, si = b + I, Sj = I + 2/(r + 1) {j = 2,...,N) 

(cf. (j3.7p ) in the integrand of ()2.ip . The integrand can then be multiphed by a suitable function 
of {ttp} to obtain an eigenvalue PDF of the form (jl.l2p . 

n 

j=2 l<j<k<n l<j<k<r{n-l) 

r{n-l) n 

^ n IIl^J -"pI'""^ '^^^2/(r+l),(r+l)iv((l-x)^). 
j=l p=l 

Integration over {Xj} corresponds to the decimation procedure D^-fi, while multiplication of 
()2.2p by this same function of {op} gives 

n 

j=2 l<j<fc<™ 

ocME2(,+i),;v((l-a;)('^+i)'+2^), 

and the sought identity results. 

Changing variables Xj i— > Xj/b in the first and second relations and taking b ^ oo gives the 
third and fourth relations. This is the Laguerre limit discussed in the sentence including (j3.14p . 
In the first identity, setting a = b, changing variables Xj ^ ^(1 — Xj)/\/2b, and then taking 
6 — > oo to obtain the Gaussian limit as discussed in the sentence including (|3.16p gives the fifth 
relation. 

The sixth relation is a consequence of Theorem [1] as it applies to (jl.9p , which tells us that 
(I2.22P is proportional to (I2.24p with Wp = e*^^ in the latter. For this set 

n = N + 1, an = b+l, aj = 1 + 2/(r + 1) (j = 1, . . . , n - 1). 

The final relation is obtained from this by setting 6 = 0, replacing by — 1, and recalling 

5 Gap and spacing probabilities 

5.1 A class of gap probabilities 

It has been revised in Section 3 that the Dixon-Anderson density (jl.ip with Sj = 1 implies 
the decimation identities for superimposed ensembles (j3.10p . (j3.1ip . These results rely on the 
eigenvalue PDF of 0E„(/) U OE„+i(/) being proportional to (l3T2]) and 0E„(/) U 0E„(/) being 
proportional to (j3.13p . No generalization of these latter facts for ensembles of the form (|1.12p is 
known. Instead let us make note of a further interpretation of the Dixon- Anderson integral in 
the case Sj = 1, or more explicitly in the case of the parameters (|3.4p with the final condition 
replaced by Sj = 1 (j = 2, . . . , A^ + 1), which does permit a generalization. 
In the latter circumstance, one has that 

„ N+l 

/ dAi • • • dXN+i n Aj (1 - Aj)' n - ^'^) 

•^^ j=l l<j<k<N+l 

= r(a + b + Ar + 2) ll"^- " H ^^^^ " "'^^ ^'''^ 

^ ' > j=l l<j<k<N 
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where A denotes the region 



1 > Ai > ai > A2 > 02 > • • • > Aat+i > 0. (5.2) 



Let 1/Cm+i be the normahzation such that the integrand of (jS.ip corresponds to the eigenvalue 
PDF 0E7v+i(a;"(l — x)*). Then with A denoting the inequahties (|5.2p we see that (jS.ip has the 
interpretation as the evaluation of the probability that there is exactly one eigenvalue in each 
interval [aj+i, a^] {j = 0, . . . , N; ao = 1, cat+i = 1), giving the result 



N 



+ ^ ^ ' ^ j=l l<j<k<N 

Similarly, replacing the condition sj = 2 in (j3.4p by sj = 1 we can obtain an analogous formula re- 
lating to the probability of the inequalities (jl.7p holding in the ensemble ME2/(^r+i),{r+i)N+r{x"'i^~ 
x)'^). However such formulas appear to be of limited interest, so we refrain from pursuing this 
further. 



5.2 Spacing probabilities 

The inter-relations of Theorem [2] have some immediate consequences for the PDF of ordered 
eigenvalues in the corresponding ensembles. To state these, let |?™'^^(A;; s; ME^ tv) denote the 
PDF for the event that in the matrix ensemble ME^ tv there is an eigenvalue at s, and exactly k 
eigenvalues to the right of s. Thus this is the PDF for the {k + l)-st eigenvalue as labelled from 
the largest. Similarly, let s; ME^^jv) denote the PDF for the event that in the matrix 

ensemble ME^ there is an eigenvalue at s, and exactly k eigenvalues between and s. In the 
case that the support of the spectrum is restricted to x > 0, this corresponds to the PDF for 
the {k + l)-st eigenvalue as labelled from the smallest. We then read off from the fourth, third 
and sixth relations in (11.161) the following result (analogous results hold for the other relations; 
however these three are representative of all situations). 



Corollary 1. One has 

p'^^"((r + l)k + r; s; ME2/(^+i),(^+i)7v+,(e-^)) 
p"^"^((r + l)k + r; s; ME^/^r+iUr+m+Ax'^e-')) 
p--(((r + l)/c + r;s;CE^/(,+i)^(,_,i)^_,, 



= p"'"^(A;;s;ME2(,+i),^(e-(^+i)-)) 

= p^'^ik; s- ME2(,+i),^(x('^+i)'^+2^e-('-+i)^)) 

= p-''^(fc;.;CE(^+_;l5+^^0 (5.3) 



Recalling (|3.25p we see that with 6 = 2/(r-|-l) the final equation in (j5.3p implies 

+ l)k + r; CE^/^^^,^^^^^,^^) = f"{k; s; CEI^,^,^^^), (5.4) 

where p^'^^"^^{k; s; CE^ tv) is the PDF for the spacing between eigenvalues which are (A; -|- l)-st 
neighbours. 

The three situations that the results of Corollary [1] are representative of the soft edge, the 
hard edge, and a spectrum singularity in the bulk. The soft edge is the neighbourhood of the 
largest eigenvalue, so called because the eigenvalue density has support beyond this region. The 
hard edge is the neighbourhood of the smallest eigenvalue, in the situation that the eigenvalue 
support is strictly zero for x < 0. The bulk is the portion of the spectrum a macroscopic distance 
from the edges (all circular ensembles satisfy this requirement), while a spectrum singularity 
corresponds to a factor of the form \x\" (for 2; ^ 0) in the one body weight. 
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Each permits a scaling in which the origin is shifted to the appropriate neighbourhood (this 
must be done in the case of the soft edge only), and the eigenvalues are scaled so that the spacing 
between eigenvalues is of order unity. For the N ^ oo limit of the ensemble ME/3^Ar(e~'^^), the 
soft edge scaling is given by [10] 

X ^ —{AN + 2{2N)^/^S(3x), (5.5) 

where > is an arbitrary length scale. The PDF for the {k + l)-st largest eigenvalue is then 
given by 

lim ^(27V)^/3g^^max(^. g. ME^,^(e-^^)) =: pf\k- s). (5.6) 

The hard edge scaling of ME^ jv(3;'^e~'^^) is given by 

P X 



2c ANsp 



(5.7) 



where s/j > is an arbitrary length scale, and the PDF for the {k + l)-st smallest eigenvalue at 
the hard edge is thus 

lim |-^p--(A;;s;ME^,^(x'^e— )) =: p^^'^k; s; a). (5.8) 
Finally, in the ensemble CE^^ the mean spacing between eigenvalues is 27r/N, and we have 

lim ^p^'^k; 27rs/N; CE^,^) =: ^^'^^'^'^•^■(A:; s; b) (5.9) 

for the PDF of the {k + l)-st eigenvalue to the right of a spectrum singularity in the bulk with 
unit density. Taking these limits in the results (15. 3p gives inter-relations between the scaled 
PDFs with P = 2/{r + 1) and p = 2{r + 1). 

Proposition 2. Choose the scale in 15. 5\) such that S2/(r+i) = (f + l)^^'^'S2(r+i); ^^'^ choose the 
scale in |5. 7\ ) such that S2/(r+i)(?' + 1)^ = ^2(^+1). One has 

Pl%.+,){{r + l)k + r;s) = p^2(r+i) (^5 ^) 
P^2/if+,){ir + l)k + r;s;a) = p^^l^^.^ik; s; (r + l)a + 2r) 
{r + l)p'^';^:lt)ii^ + l)k + r;ir + l)s;b) = p^-^^^^^- {k; s; {r + l)b + 2r). (5.10) 



Let p^^^'^^'{k; s) denote the PDF for a spacing of size s between eigenvalues which are (/c+1)- 
st neighbours, in the bulk of the circular ensemble CEjS^N scaled so that the eigenvalue density 
is unity in the N ^ 00 limit. This corresponds to the third relation in (j5.10p with b = 2/(r + 1) 
(recall the discussion leading to (j5.4p l and so we have 

+ i>2/tt?)((^ + 1)^ + + 1)^) = ^''T^^' (^-ii) 

As a consistency check on (|5.1ip it can be shown to be compatible with the asymptotic form 
(|1.17p . First one recalls that in general p^^^^'^^' {k; s) is related to {E^^^^{n; s)}n=o,i,...,k according 
to the formula (see e.g. P, Ch. 6]) 

^^"'•(^5 .) = ^ ^{k - J + l)Ef'Hr, s). 

j=0 
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For large s the term j = k will dominate, and so (jS.lip requires that for large s 

(r + l)i?2%k^i)((r + l)fc;(r + l).) ~ E'^^^l.^ik; s). 

Taking logarithms of both sides and comparing with (jl.l7p gives precise agreement with the 
latter. 
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